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Abstract 

We propose a simple, easy to implement, variant of the EXCITED method 
for variational many-body calculations for excited states. We apply this 
method to the Hybrid Multideterminant method(HMD). We test this method 
with relatively few Slater determinants by comparing the results with exact 
shell model calculations for ^^Ni using the fpd6 interaction. We obtain 
very good agreement with the exact results. 
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1 Introduction. 



In performing variational calculations for ground-state energies for many-body 
nuclear Hamiltonians, quite often we are left with the problem of obtaining rea- 
sonable values for the excitation energies. In the Hybrid Multideterminant method 
(HMD) (ref.[l]) for example, we assume that the ground-state wave function is 
well described by a sum of (usually projected to good angular momentum and 
parity) generic (that is no restrictions are imposed) Slater determinants. Both 
the Slater determinants and the coefficients of the linear combination are deter- 
mined by minimizing the expectation values of the Hamiltonian. By increasing 
the number of the Slater determinants we obtain increasingly accurate values of 
the ground-state energy as compared to the exact ones. Equivalently, after the 
minimization procedure has been carried out, we can diagonalize the many-body 
Hamiltonian in the basis of these Slater determinants and obtain a set of eigenval- 
ues. However all eigenvalues, except the ground-state, give a rather poor descrip- 
tion of the excitation energies, unless the basis consists of a very large number 
of Slater determinants. This is not very surprising since the basis has been deter- 
mined so as to minimize the ground- state energy, and therefore all other eigen- 
values are usually much higher than the exact ones. More than 20 years ago, a 
method to remove this deficiency has been proposed in variational calculations us- 
ing quasi-particle determinants, that is, the EXCITED-MAD-VAMPIR approach 
(ref.[2]). The method consists in first to obtain the ground-state wave function, 
with the appropriate number of particles, angular momentum and parity, and then 
to determine the first excited state imposing at every step orthogonality (with the 
Gram-Schmidt method) to the ground-state wave function. The method is iterated 
until all the desired excited state are obtained. We have applied this EXCITED 
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method to the HMD method. Our version of the computer program increases the 
corresponding ground-state by about 170 lines in FORTRAN 77. We found how- 
ever a much simpler method to obtain excited states which can be implemented by 
the addition of very few lines to the computer program for the ground-state. More- 
over if, after the calculation has been finished, we wish to increase the accuracy 
of the low energy wave functions, in the EXCITED method we would presum- 
ably have to repeat the calculation. The method we propose does not have this 
inconvenience and it is very flexible and easy to implement. The method is based 
on the idea that in order to increase the sensitivity of the minimization method 
on the position of the excited states energies(rather than just on the ground-state 
energy) the quantity to be minimized is some weighed average of the the first few 
low-lying eigenvalues obtained from the diagonalization of the Hamiltonian in the 
basis of the Slater determinants. For example for three states the functional to be 
minimized is F = {wiEi + W2E2 + w^E2)/{wi -1-^2 + W3) with the weights 
Wi > 0. We call this method the Centroid HMD method (CHMD). It takes very 
few lines to modify the ground- state computer program to implement this idea. 

This paper is organized as follows. In section 2 we review in more detail the 
EXCITED method (see however the original work of ref.[2]) and the centroid 
method. In section 3 we test the method by comparing the exact energies of ^^Ni 
recently obtained by shell model diagonalization in refs. [3], [4]. 



2 The EXCITED and the Centroid method. 



In the EXCITED method, that we apply here to Slater determinants rather than to 
quasi-particle states, we assume to have determined an accurate approximation to 
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the ground-state of the many-body Hamiltonian H as 



IV'(l) >=|0(1) >=E^a,l|C/a,l>, 



(1) 



where the | > are variational Slater determinants with good quantum numbers 
(angular momentum, etc.) restored with projectors. A^i is the number of Slater 
determinants necessary for good convergence to the ground-state energy (for non- 
zero angular momentum it includes also the K quantum number). The coefficients 
Qa^i are determined by solving the generalized eigenvalue problem 



El being the lowest eigenvalue. All other eigenvalues obtained in eq.(2) are dis- 
carded since they are a poor approximation to the energies of the excited states. 
The Slater determinants in eq.(l) can be generated with a two steps method. In 
the first step (addition step) Slater determinants are added to the basis one by one 
and each Slater determinant is varied so as to minimize the ground-state energy. 
In the second step (refinement step) once the basis is constructed, each Slater de- 
terminant is varied anew (ref.[5]). The motivation for this second step is that the 
first Slater determinants have been varied when the basis contained few Slater de- 
terminants. Ideally, we would want to minimize the energy by varying all Slater 
determinants simultaneously, while in the addition step only the last added one is 
varied. The refinement step corrects this partial variation. 

For the excited states one makes use of a variational ansatzs similar to eq.(l) 




(2) 



\(l^{n)>=^ga,n\Uoc,n> = 2, .., iV„), 



(3) 



a=l 



and the excited states are written as 



n— 1 



\'ll){n) >= I5nn\(i>n > + /^in|0(j) >, 



(4) 
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the coefficients /3ij — 1, ..,n) are deterniined so that < V'ilV'j >= The 
wave functions in eq. (3) are determined by minimizing 

^ < xlj{n)\H\ij{n) > 

<V'H|V'H> ■ ^ ' 

In the above equation only the ga,n\Ua.ri > are varied, all previous ga,i\Ua,i > for 
i = 1, .., n — 1 are kept fixed and the coefficients Pin (i = 1, .., n — 1) are de- 
termined by the orthogonality condition. The result of the Gram-Schmidt method 
can be recast in a compact form as follows. As in ref. [2], if we define the overlap 
matrix between previously determined non-orthogonal wave functions as 

=<0(O|0(j) > (i,i,= l,...,n-l), (6) 

and its inverse B = of dimension n — 1 and the Gram-Schmidt projector as 

n-l 

J2 <l>{i) > Bij < <f>{j)\, (7) 
then the n*'* state given by the Gram-Schmidt procedure is 

|V'(n)>=(l-^)|0(n)>, (8) 

and the quantity to be minimized is 

^ _ <0(n)|(l-^)^(l-^)|0(n)> 
< <P{n)\{l - Smn) > ' 

that is 

J-, Hnn HniBijAjn An^iBijHj^n ~\~ AniBijHjiBjj^Aj^n /-in\ 

" " A -A B -A- ■ ^ ^ 

In the above equation the sum over repeated indices, except n, is understood, and 
Hi^j —< (f){i)\H\(f){j) >. The quantities Hik — BijHjiBjk can be stored in the 
previous calculations. All overlaps and matrix elements in eq.(lO) are angular 



momentum projected, and their evaluation is certainly the most expensive part of 
the calculation. The minimization step is performed with quasi-newtonian meth- 
ods (cf. ref.[6]-[8] and references in there) once the derivatives of En with respect 
to the entries of the Slater determinants have been determined. Clearly after the 
calculation of several levels has been performed if we wish to have more accurate 
wave-functions we would have to determine anew at least some of the eigenstates. 
In the CHMD method instead we use only the ansatz 

N 

\<l>>^^ga\Ua>, (11) 

and, after the matrix elements ha',a =< U{a'\H\Ua > and Oa',a =< U{a'\Ua > 
have been evaluated we solve the generalized eigenvalue problem, written in a 
matrix notation as 

hg = ogE, (12) 

and obtain all g{a, k) for a,k — 1, .., N and the energies E^. The input to the 
minimization procedure is then F — Y^k'^kEk/ Y^Wk with positive weights Wk. 
We typically take only few states in the sum. As a rule of the thumb wi> W2> 
etc. especially for a small number of Slater determinants. The reason is that the 
location of the excited states is very sensitive to the values of ^2, ^^3, .. If we 
take comparable values for the weights, we could have an increase (rather than 
a decrease) of the ground-state energy, even though the centroid would decrease, 
because of the large contribution of the higher eigenvalues of eq.(12). Therefore 
we prefer to assign a large weight to the ground-state so a decrease in the centroid 
generally gives a decrease of all energies included in the sum. Only when the 
energies are sufficiently refined we can take comparable values for the weights. 

There is a very interesting limit of this method, namely when all Wk — 
except one. In this case we first determine the ground-state with wi — 1 (using 
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A^i Slater determinants) then we keep adding N2 Slater determinants to the basis 
setting u?! = 0,^2 = 1 (the dimension of eq. (12) is now Ni + N2), then we 
proceed setting wi — W2 — Q,W2, — 1 and add A^3 additional Slater determinants. 
Of course we do not vary the Slater determinants which belong to eigenstates 
different from the desired one. In this way we simulate the EXCITED method but 
we can improve at wish previously determined eigenstates without any concern for 
orthogonality which is automatically taken care by the diagonalization of eq.(12). 
Note that if we generate eigenstates in sequence by adding Slater determinants one 
after the other, that is if we use only the addition step without the refinement step, 
we need no new computer program at all. In fact, after A^^i Slater determinants 
have been generated for the ground-state, we can simply select the first excited 
state in eq. (12), the label of the desired eigenstate being an input parameter and 
keep adding more Slater determinants to the basis, minimizing E2 without any 
modification to the computer program. 

Note that the number of projected matrix elements to be computed in this 
limit, is the same as in the EXCITED method. When we are varying the Slater 
determinants for the excited states, the ground-state energy decreases although by 
small amounts. The main difference between this limiting case (that is all ty's set 
to except one) and the EXCITED method is that we make use of eq.(12) in a 
basis consisting of all Ni-\-N2 + N^ + .. Slater determinants, while the EXCITED 
method uses only N2 Slater determinants for the first excited state, A^3 for the 
second excited state and so on. In the next section we apply the CHMD method 
to the case of ^^Ni. 
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3 ^^NL 



As an application of the CHMD method we consider ^^Ni with an inert ^°Ca 
core, the Hamiltonian for the valence nucleons is the fpd6 Hamiltonian (ref. [9]). 
This case has been investigated recently with shell model diagonalization in the 
full space in ref. [3]. The energies for the 0i",02 and O3 states are E{Qi) — 
-203.198MeV, E{Ot) = -198.723MeV, E{Ot) = -198.204MeV (ref.[4]). For 
the 2^^ state the shell model value is E{2t) = -200.190MeV. Since we have 
only three states from the shell model calculation we consider a centroid with 
only the first three states. Let us consider first the limiting case where only one 
of the w's is non-zero. We consider the following sequence. First we generate 
15 Slater determinants with wi = 1,W2 = W3 = 0, and obtained E{Oi) = 
-203.154MeV, E{Qt) = -197.632MeV, E{Qt) = -193.330MeV. As discussed 
in the introduction, although the ground-state energy is well approximated (only 
44KeV above the exact value), the first excited 0+ state is about 1.1 MeV above 
the exact value and even worse the second excited 0"*" state is about 4.9 MeV above 
the exact value. 

We then added 15 more Slater determinants to the CHMD basis with wi — 0,W2 — 
l,ws = and obtained E(O^) = -203.157MeV, ^(0^) = -198.597MeV, 
E{Ot) = -197.881MeV. The first excited 0+ state is now only 126KeV above 
the exact value and also the energy of the O3 state has dramatically improved. We 
added then 15 more Slater determinants to the CHMD basis with wi — W2 — 
0,W3 = 1 and obtained E{Ot) = -203.160MeV, E{Ot) = -198.603MeV, 
E{0^) = -198.147MeV. The discrepancy for the 0^ state is now only 57 KeV, 
and also the energies for the Oi" and O2 states have improved. With only 45 Slater 
determinants we obtained very good values of the energies for 3 states. We de- 
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Figure 1: Energies of 0"*" states as a function of the number of Slater determinants. 

cided then to improve on the ground-state and on the first excited 0+ state. The full 
sequence of calculation is now 15 Slater determinants for 0]*" , 15 Slater determi- 
nants for O2 , 15 Slater determinants for O3", 5 more Slater determinants for the Of 
state and 5 more Slater determinants for the 0^ state. The energies thus obtained 
are E{Qi) = -203.166MeV, E{Qt) = -198.630MeV, E{Qt) = -198.132MeV. 
The discrepancies from the exact values are 32KeV,93KeV and 72KeV for the 
0^, 0^, O3" states respectively. In fig. 1 we show the behavior of the energies as 
a function of of the number of Slater determinants. The effect of the changes in 
the weights is clearly recognizable. For this set of states we considered also a 
few more calculations with all non-zero wi,W2,ws. Initially the ground-state is 
strongly overweighted, then as the number of Slater determinants is increased we 
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allow for more balanced values and at the end wi — 3,W2 — 2.5, ^3 — 1. For 
a set of 50 Slater determinants we obtained E(O^) = -203.167MeV, E{0^) = 
-198.637MeV, E{Ot) = -198.130MeV, in line with the previous results. 

For the 2+ states we evaluated the first 3 states, as before. We have used 50 
Slater determinants and used non-zero values for the weights. Our results are 
E{2t) = -200.122MeV, E{2t) = -198.064MeV and E(2j) = -196.725MeV. 
The energy of the 2^ state is about 70KeV higher than the shell model value. The 
shell model values for the other states are not available. Although, there is still 
some room for a further decrease in the energies, the excitation energies seem to 
be stable enough as a function of the number of Slater determinants, whithin an 
uncertainty of few tens of KeV. 

In conclusion, we have presented a simple method to evaluate the non-yrast 
eigenstates with minimal, if any, modifications to the ground- state computer pro- 
grams. 
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